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We apply a recently developed functional renormalization group (fRG) scheme for quantum spin 
systems to the spin-1/2 antiferromagnetic XXZ model on a two-dimensional square lattice. Based 
on an auxiliary fermion representation we derive flow equations which allow a resummation of the 
perturbation series in the spin-spin interactions. Spin susceptibilities are calculated for different 
values of the anisotropy parameter. The phase transition between planar and axial ordering at 
the isotropic point is reproduced correctly. The results for the critical scales from the fRG as 
quantitative measures for the ordering temperatures are in good agreement with the exact solution 
in the Ising limit. On the easy-plane side, the deviations from critical temperatures obtained with 
quantum Monte Carlo are larger but still acceptable. However, at the isotropic point the Mermin- 
Wagner theorem is violated such that a precise description of the behavior in the vicinity of the 
phase transition is not possible. We discuss possible reasons for these discrepanies. 
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I. INTRODUCTION 

Renormalization group (RG) methods for fermions de- 
veloped over the last decades have become a widely used 
tool for correlated electron systems^^. Besides a quali- 
tative understanding of fundamental low-energy aspects 
of weakly interacting many-fermion systems such as the 
stability of Fermi liquidsSEI or the effective behavior of 
impurities in Luttinger liquids^!, RG methods have also 
been used extensively for investigating competing order- 
ing tendencies in two-dimensional lattice systems or the 
different r egim es in impurity problems in one or zero 
dimensions^'. Functional RG (fRG) methods do not 
only concentrate on a smaller number of possibly rele- 
vant terms in the theory, but aim to resolve the flow of 
the generating functional of the theory or the respective 
vertex functions in as much detail as possible. This al- 
lows one to monitor how effects at intermediate energy 
scales drive the flow to certain fixed points, and at which 
scales typical low-energy behaviors such as power-laws or 
scaling should actually hold. While for bosonic problems 
one can implemen t th e fRG non-pert urbativel}^^, for 
fermionic problem d 4 ' 10 ! one has to work with a truncation 
in the powers of fermionic fields. This usually restricts 
the validity to weak to moderate interaction strengths. 
Many applications use the fRG method primarily as a 
qualitative tool for the exploration of theoretical path- 
ways to different low-energy behaviors. Nevertheless, it 
is important to note that for most cases alluded to above, 
the approach is controllable in the limit of small interac- 
tion strengths. For example, the exponents in Luttinger 
liquids obtained using the fRG are consistent with non- 
perturbative or exact solutions^. Likewise, the <i-wave su- 
perconducting pairing instability in the two-dimensional 
Hubbard model can be seen to be immune against cor- 



rections in the limit of small interactions^. 

Competing orders and strong quantum fluctuations are 
also encountered in quantum spin systems^ESl. I n this 
very active research area, analytical or field theoretical 
approaches are however even more involved than for mod- 
erately interacting fermions due to the constraint of a 
fixed total spin per lattice site. Direct numerical meth- 
ods are often the most reliable approach to these strongly 
correlated problems, but finite-size effects and sign prob- 
lems due to frustration effects can make it difficult to 
obtain clear results. Hence new methods might help. Re- 
cently, the fRG has been formulated for auxiliary 5=1/2 
fermions describing localized quantum spins. In a series 
of papers^N22'it was shown that this approach gives qual- 
itatively appealing results e.g. for frustrated systems on 
various two-dimensional lattices. As known from other 
problems, one advantage of the fRG with respect to other 
numerical methods is the flexibility in treating different 
lattice geometries and dimensions. With manageable nu- 
merical effort one can study rather large, frustrated sys- 
tems without encountering sign problems. The formula- 
tion is transparent; for example, competing trends can be 
identified and manipulated (by modifying terms) in the 
calculations. Hence, it is an interesting question how far 
the spin-fRG approach can be pushed and how it com- 
pares quantitatively to other methods, in cases where 
benchmarking is possible. 

Testing the spin-fRG is also interesting from a theo- 
retical point of view, for several reasons. Expressing the 
localized spins in terms of auxiliary fermions brings some 
important differences with respect to 'normal' fermions 
on a lattice, (i) As mentioned above, the auxiliary 
fermion number on each site should remain constrained 
to a fixed value in order to guarantee a faithful represen- 
tation of the localized quantum spins. Most of the pre- 



vious applications of the spin-fRG used a mean-field-like 
approximation at T — where the constraint is treated 
on average. We will adopt this approximation here as 
well. In Ref. RBI the authors also report on the possibility 
to go beyond this approximation by using an imaginary 
chemical potential at nonzero temperature T. How- 
ever, this implies a substantial increase in the numeri- 
cal effort and we do not consider this option here, (ii) 
The bare action in terms of the auxiliary fermions does 
not contain a kinetic energy, i.e. the auxiliary fermions 
are dispersionless. This is, of course, a prerequisite for 
the particle- number constraint. Fortunately, in the fRG 
for this system the dispersionless character is preserved 
and the auxiliary fermion self-energy remains local, which 
also significantly simplifies the numerical treatment. On 
the other hand, the absence of a kinetic energy means 
that, at least at T — > 0, the system is strongly cou- 
pled, and a priori there is no reason why a perturbative 
fermionic fRG setup (or any other theory perturbative 
in the interactions) should work well, (iii) Quadratic in- 
teractions between localized spins translate into quartic 
auxiliary fermion interactions that are bi-local in the lat- 
tice index, i.e. consist of two pairs of fermion fields with 
the same site index. Again, this structure is preserved 
under the fRG flow. This is a tremendous simplification 
with respect to the 'normal' fermionic case where in gen- 
eral also interactions with more than two site indices oc- 
cur. Concluding this discussion of the auxiliary fermionic 
description one can say that points (i) and (ii) are at 
least two major uncertainties in the spin-fRG approach. 
On the positive side, the simplified structure of both the 
bare and the flowing theory allows one to capture basi- 
cally the full frequency dependence of the self-energy and 
the four-point interaction vertices, and in addition makes 
it possible to include the self-energy feedback in the flow 
equations without too much effort. This means that in 
the spin-fRG, one reaches an advanced approximation 
level (for extended systems in more than one dimension) 
that has before o nly be en used for impurity systems in 
very selected case d 22 * 2 ^ . Hence, it is a very interesting 
question how this combination of positive and negative 
expectations performs if one tries to compare the results 
with other techniques. 

The goal of the present paper is to provide a quanti- 
tative test of the spin-fRG. In order to provide a well- 
understood testing ground, we study the spin-1/2 an- 
tifcrromagnetic XXZ-model, which as function of the 
anisotropy parameter is known to exhibit two differ- 
ent orderings, and where the critical temperatures for 
these orderings show a nontrivial variation. The critical 
scales A c for these orderings computed with the spin- 
fRG are compared with the critical temperatures T c ob- 
tained from quantum Monte Carlo (QMC) calculations. 
We note that while it is in principle possible to imple- 
ment the spin-fRG for finite temperatures and hence to 
directly measure T c , this would drastically increase the 
numerical effort. Thus in practice one aims to relate the 
energy scales A c found more easily at T = to the ac- 



tual T c . Our work provides insights to where and to what 
extent this is possible. 

This paper is organized as follows: after introducing 
the anisotropic Heisenberg model, the auxiliary-fermion 
representation and the technical implementation of the 
fRG specific to spin systems is shortly reviewed. We 
then present the results for the spin susceptibility and 
the critical scales of the fRG as estimates for the critical 
temperature and discuss the merits and limitations of the 
spin-fRG. 



II. MODEL 

In the following, we consider the spin-1/2 XXZ model 
with antiferromagnetic nearest-neighbor exchange inter- 
action J > 0, defined by the Hamiltonian 

// ./^(V'V ■ S»S? . A.S>;). (1) 

Here, Si denotes the spin-1/2 operators on lattice site i of 
a two-dimensional square lattice with periodic boundary 
conditions. Furthermore, A is the exchange anisotropy, 
which we restrict to the antiferromagnetic domain, i.e. 
A > throughout this study. For the above model, both 
the ground state properties as well as the finite tempera- 
ture phase diagram are well established, cf. e.g. Ref. l24l 
Here, we mention those features that are relevant for our 
assessment of the spin-fRG method. 

In the easy-axis region, A > 1, the ground state is 
Neel-ordered along the longitudinal (z) direction, while 
within the easy-plane region for A < 1, the ground state 
exhibits long-range transverse antiferromagnetic ordering 
within the xy-plane. Within the easy-axis region, long- 
ranged longitudinal Neel-order emerges already below a 
finite transition temperature T c > 0. This transition 
corresponds to the breaking of the discrete Zi symme- 
try of H, and thus belongs to the universality class of 
the two-dimensional Ising model. In particular, in the 
large- A limit, A —¥ oo, the transition temperature scales 
like T c 2.269A/4, as obtained from the exact solu- 
tion of the two-dimensional Ising modeP^ (the factor of 
1/4 results from the different normalization of the spin 
variables usually employed in the classical Ising model 
as compared to the spin-1/2 operators entering H). The 
ground state ordering in the easy-plane region, i.e. for 
A < 1, relates to the breaking of the U(l) symmetry 
of H, such that in two dimensions the transverse long- 
range order does not persist to nonzero t emp eratures, 
consistent with the Mermin- Wagner theoreirP2. Still, the 
system e xhibit s a Berezinskii-Kosterlitz-Thouless (BKT) 
transitio n 1 27 ! 28 ! at a critical temperature Tbkt > driven 
by the unbinding of vortex excitations above Tbkt out of 
a low-temperature quasi-long-ranged ordered phase with 
algebraic decay of the transverse spin correlation func- 
tion. This phase is furthermore characterized by a fi- 
nite spin stiffness, which exhibits an universal jump at 
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Tbkt- At A = 1, the Hamiltonian H exhibits the full 
SU(2) symmetry and the antiferromagnetic order is con- 
strained to zero temperatures, such that at this isotropic 
point the critical temperature vanishes, i.e. T c = 0. The 
A-dependence of both the Ising transition temperature 
and the BKT transition temperature in the vicinity of 
A = 1 has been analyzed in Ref. |2T)1 

However, in order to assess the quality of the spin-fRG 
approach to the XXZ model, we require an estimation 
of the critical temperature T c for a wider range of differ- 
ent A-values. Since the model is not frustrated, we can 
indeed employ numerically exact, unbiased QMC meth- 
ods for this purpose. Here, the critical temperatures in 
both the easy-axis and easy-plane regions have been ex- 
tracted from QMC simulations using the stochastic se- 
ries ex pansio n method^ with generalized directed loop 
updates^ 1221 . In particular, we simulated finite lattices 
with N = L 2 spins and with linear system sizes L rang- 
ing up to 32 using periodic boundary conditions. From 
the finite size data, the critical temperature for the Ising 
transition was obtained from a finite-size analysis based 
on the critical Binder rations. Within the easy-axis re- 
gion, the spin stiffness was calculated from measuring the 
spin winding number fluctuations^. From the universal 
stiffness jump at the BKT transition, the transition tem- 
perature was then obtained^. The values of the tran- 
sition temperatures for |A — 1| < 0.1 can also be taken 
from Ref. |29] and for A = from Ref. US The results 
of these calculations will be presented below. However, 
before presenting these data, we provide a short intro- 
duction to the spin-fRG method in the following section. 



III. AUXILIARY FERMIONS AND SPIN-FRG 

Here we briefly describe the spin-fRG scheme intro- 
duced by Reuther and Wolfle^. In order to be able to 
implement the fermionic fRG scheme^^ for the quan- 
tum spin system, one chooses a representation of the spin 
operators in terms of Abrikosov auxiliary-fermiona^. 
This representation requires two fermionic operators /, Q , 
a =t)4o for each lattice site i, 

St = \Y,fl<^ ( 2 ) 

a/3 

with tr M being the Pauli matrices. By construction, the 
representation satisfies the correct commutation relation 
of spin operators. The advantage of this representation 
lies in its quadratic form which allows for the applica- 
tion of Feynman-diagram techniques. However, the in- 
troduction of auxiliary fermions implies an enlargement 
of the Hilbert space, as a single site can now also carry 
total spin zero, either by being empty or by being dou- 
bly occupied. In order to obtain a faithful representation 
of the quantum spin-l/2-system, these unphysical states 
have to be projected out. This is formally achieved by 



the requirement that the single-occupancy operators 

Qi = (3) 

a 

all possess the eigenvalue 1 on any physical state |phys) 
occurring in the thermodynamic averages of the fermionic 
system, i.e. Qi|phys) = |phys). Thus all unphysical 
states are discarded. 

The single-occupancy constraint immensely compli- 
cates the fRG treatment. A convenient approximation 
is to enforce the constraint only on average by requiring 
(Qi) — 1- With this, of course, local particle number 
fluctuations are still possible. Due to the particle-hole 
symmetry in the dispersionless pseudo-fermion system 
(Qi) = 1 is fulfilled for /i — 0. For an exact projec- 
tion without fluctuations, Popov and Fedotov^ proved 
the mutual cancellation of the unphysical contributions 
at each lattice site in the thermodynamic averages if an 
imaginary chemical potential fi — —inT/2 is used. How- 
ever, as the limit T —¥ and the suppression of fluc- 
tuations at T — do not necessarily commute, the two 
procedures need not to be equivalent at T = 0. Following 
Ref. [2]we nere consider the average projection (Qi) = 1 
and vanishing temperature T = for simplicity. More 
details to this point can be found in Ref. |3U 

After transforming the initial spin-system into a 
fermionic one, the requirements to apply the fermionic 
fRG are fulfilled (except for the lack of a small expansion 
parameter, see the discussion below). The ideePis to sum 
up the contributions to the 1-particle irreducible (1PI) 
vertex functions from high to low energies. To achieve 
this, we introduce a cutoff A into the bare Green's func- 
tion which is replaced by 

g o,a m = e(H-A) 

Inserting G°' A into the generating functional of the 1PI 
vertices and performing the derivative with respect to A 
leads to an infinite hierarchy of flow equations. For nu- 
merical calculations this hierarchy has to be truncated. 
It turns out that a truncation after the four-point vertex 
in the straightforward form of the 1PI fRG hierarchy^! 
is not sufficient to adequately describe the competition 
between order and disorder fluctuations^. An improved 
scheme suggested by Katanin 40 inclu des m ore self-energy 
corrections such that Ward identifie d 41 * 42 ! are better ful- 
filled. With this Katanin-modification, which we will use 
as well, the results for the J\- J2-model were founds to 
be in good agreement with other studies. However, we 
will see that the spin fRG still overestimates the mag- 
netic ordering tendencies during the flow. Unfortunately, 
due to the increasing numerical effort the inclusion of 
higher-order contributions beyond the Katanin trunca- 
tion scheme is not feasible. 

Within the Katanin-modified lPI-scheme, the RG 
equations for the self-energy S A and the 2-particle vertex 
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r A read 
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— r A , 
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(6) 



where the multi-indices 1 = W\,i\,a\ include the fre- 
quency, site and spin, and the sum over 1 contains an 
integral over u>\ as well as a sum over i\ and ol\. Fur- 
thermore, 



is the full propagator and 



S A (ioj) 



e(M-A) 

icu — T, A (iuj) 



iw — E A (zw) 



(7) 



(8) 



denotes the single-scale propagator. The self-energy is a 
purely imaginary odd function of the frequency and can 
be written in the form 



e a m = -z 7 A ( w ). 



(9) 



We observe that in the diagrammatic expansion the 
Green's function remains strictly local— and the momen- 
tum dependence of the susceptibility is generated by the 
non-local exchange couplings. If one uses the symmetry 
relations of the system and energy conservation, one can 
parametrize the 2-particle vertex in the form 

+ T t y i 1 i 2 (^^ u )« 1 ,a 1 < 2 ,a 2 + (r l 1 ,a 1 < 2 ,a 2 ) 
~^~^di 1 i 2 ( S ' ^' !l )^a 1 ;ai^a 2 /a2 ^i 1 /ii^i 2 '»2 
r « lJ2 ( s 1 "^) cr a 1 ,a 2 cr a 2 ,a 1 

+r A ■ (s u t)(a x a x +a y a v ) 

< 1 xyi 1 i 2 \°i ' l /ro 1 iaj u a 2 /ai ~ u a x i a 2 u a 2 i ai J 
+^di 1 i 2 ( S ' M ' t)$a 1 ,a 2 fia 2 ,a 1 x Si 1 ,i 2 Si 2 ,i 1 (10) 

with 'Mandelstam'-variables s = w\ + u>2, t = ujy — u)\ 
and u = uj\ —LU2> ■ To be more precise, the vertices depend 
only on the difference \ix — z 2 1 due to translations! invari- 
ance. The resulting RG equations for j(uj), T Z ij(s,t,u), 
T xy ij (s, t, u) , and Tdij (s, t, u) including the full frequency 
dependence are explicitly stated in Ref. |43l The initial 
conditions are 



7 ■ 



= 0, 



r A= °° (s t u) = < 4 



for |ii-i 2 | = l, 
otherwise, 



(11) 
(12) 



for \i x - i 2 \ = 1, 
otherwise, 



r A =™( s ,t,u) = o. 



(13) 
(14) 



We note that the spatial structure of the interactions only 
enters the initial conditions. 

We can also drop some terms in this 'full' set of fRG 
equations and thereby reproduce single-channel sum- 
mations of random-phase-approximation (RPA) charac- 
ter. Below we will present RPA data from such a re- 
stricted flow, where only specific terms in the i-channcl 

(~ r i\4;l,3 r 3,2';4.2 + ( 3 ° 4 )) ° f tllC fl ° W e q uation Q 

for the 2-particle vertex are taken into account, as ex- 



plained in Ref. 14] Inserting the representation (10 1 into 
the i-channel gives various terms. The terms involving 
T A , , remain zero, as the initial condition r4 = °° — 
is conserved in the flow with the f-channel only. The 
remaining terms can be analyzed most easily by Fourier- 
transforming the vertices T A , . . from i-i — in to q. In 

° z ( xy i\i 2 

reciprocal space, the vertices only depend on one wave- 
vector q. In the i-channel for the flow of T A , -. there are 
now terms that do not couple different qs, while other 
terms contain an internal sum over another wave-vector 
<f . The first class of terms can be interpreted as bubble- 
like RPA, while the second class could be called vertex 
corrections, and are dropped in the RPA-approximation. 
This leads to the RPA-like flow equation for the vertex, 
explicitly written out in Eq. (40) in Ref. Q3J We discuss 
two versions: In the "RPAO" the self-energy is neglected, 
whereas for the "RPA+" self-energy contributions are 
taken into account as in the usual self-consistent RPA 
without cutoffs. For the flow equation of the self-energy 
no further selection of the diagrams is necessary, here 
only a Fock-like contribution remains. Note that in spe- 
cific fermionic models, integrating the RPA + H artree 
terms gives results equivalent to mean- field theorjHSIM] 
but RPA + Fock (or, equivalently, mean-field+Fock, as 
the Hartree term vanishes) as in our case does not allow 
such an equivalence. Furthermore, note that the self- 
consistent RPA for the 2D Heisenberg-antiferromagnet 
containing the same class of RPA + Fock diagrams has 
been found to prevent long-rang order at T > in ac- 
cordance with the Mermin- Wagner theorem, but at the 
price of overdamping the spin-wave excitations at small 
frequencies^. When comparing this to the RPA+, it is 
important to note that while diagrams with the same 
topology are used, the propagators on the internal lines 
differ, and in general the fRG results integrated down 
to zero scale cannot be expected to match those of self- 
consistent schemes. 



IV. RESULTS 

We have solved the RG equations ^ and ^ with 
the initial conditions (11|-(14| numerically. Following 



Ref. [T3] we used a logarithmic frequency discretization of 
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Nf = 30 positive frequencies, if it is not explicitly stated 
otherwise, with uj m i n = 0.001J and w m ax ~ 250J This 
ensures that the small frequencies are covered in detail 
and some frequencies higher than A x J are used. Neg- 
ative frequencies are implicitly included via symmetry 
relations. Unless stated otherwise, all results were ob- 
tained for a lattice with 10 x 10 sites. Different from 
the original treatment by Reuther and Wolfle^, the in- 
teractions are not truncated in range. Instead periodic 
boundary conditions are used. 

A. Susceptibility 

First we compute the static longitudinal spin-spin cor- 
relation function xfj by a diagrammatic calculation 

/>oo 

X?/(i" = 0)=/ dT(T T [S?(r)S*(0)]) (15) 
Jo 

and obtain the Neel susceptibility x zz through 

X zz [q= (n,n),uj = 0] = ^(-l)^+^ X o|(^ = 0) , 

where j x (j y ) denotes the distance in x(y)-direction to 
an arbitrarily chosen reference point. The transverse 
spin-spin correlation function xfj i s defined analogously. 
These quantities only have a strict physical interpretation 
at A = 0, but we will consider them for finite A as well 
and interpret a divergent flow as an indicator for a mag- 
netic instability of the system. Again, using the physical 
temperature as a cutoff or working at T > would make 
this interpretation more straightforward, but would also 
strongly increase the numerical effort. The experience 
from other fermionic systems^ shows that divergences 
with respect to A and T are qualitatively equivalent. 

The results for different anisotropies A are shown in 
Fig. [I] For A < 1 we find a dominating transverse 
susceptibility x xx '■> whereas for A > 1 the longitudinal 
susceptibility x zz diverges at a larger scale A. At the 
isotropic point A = 1 both susceptibilities are identical 
as required by the full SU(2) invariance. This finding 
thus reproduces the well-known^ phase transition be- 
tween planar and axial ordering at A = 1. Furthermore, 
it has to be stressed that due to the generic divergence of 
the RG flow at a finite scale A c it is impossible to obtain 
quantitative results for the susceptibility which could be 
compared to other methods on an absolute scale. 

B. Critical scales and temperatures 

In Ref . H"9l Reuther et at. directly interpret the scale A c , 
at which the fRG flow diverg critical temperature 

T c for the phase transition into the ordered phase. In 
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FIG. 1. (Color online) Flow of the static transverse (x xx ) 
and longitudinal (x zz ) susceptibility for different anisotropy 
parameters A. The results confirm the existence of a phase 
transition between planar and axial ordering at A = 1. 



order to test this interpretation we have first studied the 
dependence of A c on the system size and the frequency 
discretization. The results shown in Fig. [2] suggest that 
the scale A c is indeed rather constant at least in the range 
of lattice sizes studied here. A mild finite-size effect is 
however visible and indicates that A c tends to increase 
with system size. As one can see in Fig. [3j for a large 
enough number of frequencies Nf the scale A c is essen- 
tially independent of the discretization. In contrast to 
the system size dependence no monotonous dependence 
on Nf is observed. Still for larger A, the relative differ- 
ence in A c for different numbers of frequencies becomes 
smaller. 

As a next test we analyze the Ising limit (A — > oo) for 
which the scale A c is shown in Fig. [4] As discussed in 
Sec. |llj the critical temperature T c = 2. 269 A/4 is known 
exactly in this limit. We observe that as a function of 
large A the values of A c of the RPA0, the RPA+ and the 
spin-fRG converge, but that the limiting values differ. 
Comparing the result for the RPA0 and RPA+ approxi- 
mations we conclude that the inclusion of the self-energy 
improves the results, but that both methods still overes- 
timate the ordering tendencies of the system. The spin- 
fRG presents an improvement over the RPA as the calcu- 
lated A c agrees rather well with the critical temperature 
T c in the Ising limit. From this finding we may conclude 
that A c can indeed be interpreted as a critical tempera- 
ture A c = T c (at least in the Ising limit). However, as we 
show next, this interpretation becomes less satisfactory 
when leaving the regime of strong Ising anisotropy. 

In particular, we consider the regime A w 1 where 
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FIG. 2. (Color online) Critical flow parameter A c extracted 
from the divergence of the spin-fRG flow as a function of 
the anisotropy parameter A for different system sizes. A c 
increases with the system size. 
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FIG. 4. (Color online) Critical flow parameter A c as obtained 
from different approaches in comparison to the exact result for 
the critical temperature in the Ising limit. In comparison to 
the RPA approaches the spin-fRG reproduces the Ising limit 
far more correctly. The errors of the QMC data are smaller 
than the size of the symbols. 
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FIG. 3. (Color online) Critical scale A c for different numbers 
of positive frequencies Nf used in the discretization. 



we compare our results with the critical temperature de- 
termined from the QMC analysis. As one can see in 
Fig. [5j the spin-fRG improves the RPA in this regime as 
well, but the deviation from the QMC results is clearly 
larger than in the Ising limit. Nevertheless, except for 
the direct vicinity of the phase transition at A = 1, the 
spin-fRG scale A c correctly describes the dependence of 
T c on the anisotropy, and for A > 1.5 the fRG critical 
scales A c provides a reasonable quantitative estimate for 
T c , whereas for A < 0.5 the difference between fRG and 
QMC is significant. 

The strongest deviations between the critical scale ex- 
tracted from the spin-fRG and the QMC result are found 
close to the isotropic point A = 1, where the fRG yields 
a finite critical temperature in violation of the Mermin- 
Wagner theorenPfil. Nonetheless the A-dependence of T c 
obtained in the spin-fRG shows a kink-like feature at 
A = 1, which is not observable in the RPAO calcula- 
tion and is only very weak in the RPA+. One may ar- 



FIG. 5. (Color online) Critical flow parameter A c as a func- 
tion of the anisotropy parameter A near the isotropic point. 
The errors of the QMC data are smaller than the size of the 
symbols. 



gue that at T = ordering indeed occurs and hence the 
nonzero critical scale obtained from spin-fRG does not 
violate the Mermin- Wagner theorem. We note, however, 
that the isotropic Heisenberg model is gapless and it is 
not clear which energy scale A c thus represents. A more 
plausible way to understand the data is that the trun- 
cated spin-fRG fails to capture a significant part of the 
physics of collective fluctuations that would be needed to 
restore T c = 0. 



C. Auxiliary fermion self-energy 

Closer inspection of the RPA data for A c reveals that 
RPA+ shows a (weak) kink at A = 1 while the RPAO 
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data are monotonically increasing. The main difference 
between the two methods is the generated self-energy of 
the auxiliary fermions. Thus a more detailed analysis of 
this quantity in the spin-fRG may lead to some insights. 



The auxiliary fermion self-energy generated in the 
spin-fRG is not a physical quantity, since the auxiliary 
fermion Green's function connects physical and unphys- 
ical sectors of the fermionic Hilbert space. Neverthe- 
less, analyzing the self-energy sheds some light on the 
shortcomings of the spin-fRG and the differences to other 
methods. At the end of the flow, near A c , the frequency- 
dependent self-energy 7 A (w) shows a peak structure (see 
Fig. [6b, where the maximum occurs at uj w A at least 
for A > 1. In this regard, the fRG self-energy qualita- 
tively differs from the one found in self-consistent RPA^, 
which has a peak near zero frequency. The delayed gen- 
eration of the self-energy in the fRG may be the reason 
why the self-consistent RPA of Ref. H5| and the RPA+ 
implemented here do not agree although the same class 
of diagrams is kept in both approaches. 



For a two-site system the exact auxiliary fermion self- 
energy can be calculated^ at vanishing chemical poten- 
tial from a Lehmann representation yielding 7 ~ 1/uj. 
However, for this case the self-energy generated by the 
fRG possesses the same form as for the 2D model con- 
sidered here (see Refs. H51 and ITFj). i.e. with a peak at 
nonzero frequency whose position depends on A. This 
difference possibly originates from an incorrect descrip- 
tion at small frequencies. Based on this hypothesis, we 
can give an argument why the spin-fRG critical scales 
agree with the QMC results for T c in the Ising limit. 
Here, the flow diverges already at larger scales A. As 
the self-energy enters the flow equations via the propa- 
gators ~ + 7(k>)], its form for to > A only weakly 
affects the flow. However, for Awl, the flow reaches 
smaller scales and this influence becomes more and more 
important. This qualitatively explains the good agree- 
ment in the Ising limit and may to some extent (more 
discussion is given further below) account for the failure 
of the method to correctly describe the critical tempera- 
ture near the isotropic point and, in particular, the vio- 
lation of the Mermin- Wagner theorem. Here a stronger 
(or singular) self-energy at small frequencies would help 
to avoid the divergence of the flow at a non-zero scale. 
We note that in self-consistent studies in the auxiliary 
fermion language^ the results could be improved by an 
artificially introduced stronger suppression of the auxil- 
iary fermion spectral function at small frequencies. Sim- 
ilarly in the spin-fRG, a manipulation towards a larger 
self-energy at small frequencies would lead to a critical 
temperature T c — at the Heisenberg point and therefore 
the fulfillment of the Mermin- Wagner theorem, but the 
critical temperature especially in the easy-plane regime 
would be strongly affected and vanishing, too. 
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FIG. 6. (Color online) Generated self-energy at the end of the 
fRG flow, i.e. at A = A c , for different values of A. During the 
flow the peak builds up and moves in from higher frequencies. 



V. DISCUSSION 

The goal of this paper was a comparison of spin-fRG 
results with those of other methods to perform a quali- 
tative and quantitative assessment of the spin-fRG. We 
were able to confirm the phase transition between pla- 
nar and axial ordering at the isotropic point A = f. 
The identification of the cutoff scale A c with the criti- 
cal temperature T c yields rather accurate results in the 
Ising limit, but at the same time leads to a violation of 
the Mermin- Wagner theorem at the isotropic Heisenberg 
point and to deviations in the easy-axis and, to larger 
degree, in the easy-plane regime. The generated auxil- 
iary fermion self-energy seems too weak to prevent diver- 
gences in the fRG flow for A c > 0. Possible reasons for 
this error are: 

(i) The long-wave physics is not included as the system 
is still finite. We cannot exclude that the spin waves are 
just too constrained by our limited system size. We note, 
however, that for the available system sizes we observe 
an increasing rather than decreasing A c with system size 
(see Fig. §. 

(ii) The mapping from the spin- to a fermionic sys- 
tem, where the constraint is fulfilled only on average 
and does not apply to quantum fluctuations. Here it ap- 
pears problematic that the auxiliary fermion self-energy 
remains finite at small frequencies, in contrast with the 
exact solution of the the two-site problem, and possibly 
in contrast with the requirement of zero particle-number 
fluctuations of the auxiliary fermions. The finiteness of 
the self-energy is a consequence of the structure of the 
frequency flow, and major changes to the formalism may 
be needed to alter this behavior. This projection problem 
could also be investigated more systematically by includ- 
ing a nonzero temperature, where an exact projection 
scheme using an additional imaginary chemical potential 
is available^. The major disadvantage consists in the 
fact that the imaginary chemical potential breaks vari- 
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ous symmetries increasing the numerical effort substan- 
tially. For the Heisenberg model such calculations have 
been performed in Rcf. 39 As the results showed no 
qualitative differences to the average projection scheme, 
it seems that this is not the main source of error. 



(iii) The RG treatment of the system with a trunca- 
tion after the four-point vertex is not controlled because 
there is no small expansion parameter in the XXZ model. 
In bosonic descriptions of 0(3)-models, the four-boson 
coupling is essential to drive the critical temperature to 
zero and hence to fulfill the Mermin- Wagner theorem. 
In fermionic language, this coupling would correspond 
to a 1PI eight-fermion vertex that is not present in the 
truncation used here for the spin-fRG. In this sense, the 
interaction between low-lying spin fluctuations that sup- 
presses also the zero-frequency ordering might not be im- 
plemented correctly in this approach. This is also sup- 
ported by the observation that the critical scales deviate 
more strongly from the QMC results for T c on the easy- 
plane side where the spin waves are gapless. 



Furt her te sts have been performed for one-dimensional 
modelfPEl. Here again the spin-fRG performs better 
than the RPA. For example in the latter there still exists 
a finite scale A c at which the flow breaks down, while 
in the former the RG flow can be performed down to 
the smallest frequency uj m in- The spin-fRG correctly re- 
produces qualitative features like antiferromagnetic cor- 
relations. However, on a more quantitative level, the 
power-law decay of the spatial correlations in the Heisen- 
berg chain is not reproduced by the spin-fRG, which 
rather predicts an incorrect exponential behavior. Simi- 
lar shortcomings are seen in other models, e.g. the spin- 
fRG generically yields too short-ranged correlations in 
both gapless and for gapped systems. Furthermore the 
ground-state energies and their dependence on the sys- 
tem parameters and number of lattice sites are not repro- 
duced. Hence, the spin-fRG cannot be regarded as a reli- 
able theoretical method for the study of one-dimensional 
spin system. 



VI. CONCLUSION 

In conclusion, the spin-fRG considerably improves 
RPA results for low-dimensional quantum spin systems. 
The qualitative behavior of the spin correlations on short 
scales is reproduced correctly. For the two-dimensional 
XXZ model studied here, the spin-fRG does detect the 
quantum phase transition at A = 1 and the critical scale 
A c can be used as an estimate for ordering temperatures 
away from the isotropic point. For all anisotropics 
| A — 1| > 0.1, the deviations in A c from T c are less than 
a factor of 2. At the isotropic point A = 1, however, 
the spin-fRG predicts a finite ordering temperature in 
violation of the Mermin- Wagner theorem, which also 
spoils the quantitative validity of the spin-fRG results in 
the vicinity of A = 1. Furthermore, as discussed above, 
due to the generic divergence of the RG flow at a finite 
scale A c it is impossible to obtain quantitative results for 
the susceptibilities or correlation functions which could 
be compared to other methods on an absolute scale. 
These two drawbacks are well known from the fRG 
treatment of itinerant two-dimensional many-fermion 
systems^, and also the appealingly advanced approxi- 
mation level of the spin-fRG with the full-self-energy 
feedback and frequency-dependence cannot remedy this 
deficiency. Hence, one may conclude that the spin-fRG 
should be rather used to explore phase diagrams in 
situations where strongly differing types of ground 
states are expected to compete, i.e. ground states that 
differ in their correlations already on short distances. 
Such situations were actually addressed in most of the 
spin-fRG papers up to now. 
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